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Abstract 

We show that the tilted perfect fluid Bianchi VIo family of self-similar models found 
by Rosquist and Jantzen [K. Rosquist and R. T. Jantzen, Exact power law solutions of 
the Einstein equations, 1985 Phys. Lett. 107A 29-32] is the most general class of tilted 
self-similar models but the state parameter 7 lies in the interval f |, |j. The model has a 
four dimensional stable manifold indicating the possibility that it may be future attractor, 
at least for the subclass of tilted Bianchi VIo models satisfying = in which it belongs. 
In addition the angle of tilt is asymptotically significant at late times suggesting that for 
the above subclasses of models the tilt is asymptotically extreme. 
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In cosmological context transitively self-similar Bianchi models (spatially homogeneous models 
admitting a four dimensional group of homotheties acting simply transitively on the spacetime 
manifold) keep a privilege position among other high symmetric models (i.e. spacetime homo- 
geneous models which admit a simply transitive group of isometries or solutions with special 
geometric or algebraic properties). This is due to the fact that transitively self-similar Bianchi 
models arise as equilibrium points in the dynamical state space of more general Bianchi models 
and determine various stable and unstable submanifolds which may be identical to some of 
the invariant sets of the orbits, consequently providing a way to gain deeper insight into their 
asymptotic behavior pQ. At the same time the assumption of self-similarity has the advantage 
of reducing the Field Equations (FE) to a purely algebraic (closed) form which, if not can be 
solved explicitly, is used for doing numerical simulations in Bianchi cosmologies. 
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Up to now, the problem of finding the complete set of self-similar Bianchi perfect fluid 
models with linear equation of state of the form p = (7 — l)p, is still open due to the (highly) 
non-linear character of the resulting algebraic system. Considerable progress has been made 
for vacuum and non-tilted perfect fluid models (orthogonal spatially homogeneous cosmologies 
|2|) in which the whole class of self-similar models has been determined [3J. In the case of 
tilted perfect fluid models the situation is more complex and only some particular results are 
known. For example Bradley has stated that there do not exist tilted dust self-similar models 
whereas it was proved recently that irrotational Bianchi type VIo and models of type VIIo, VIII 
and IX are not self-similar [5J E] • 

The existent self-similar solutions have not been fully exploited apart from the case of 
Bianchi type II models in which it has been proved [7j that self-similar models may act as 
future attractors for more generic models. In the general Bianchi type II self-similar solution 
found by Hewitt jH] the state parameter lies in the interval 7 G (y, 2j. This self-similar model 
has zero vorticity and non-zero acceleration and was used to analyze in full detail the whole 
class of type II models. In fact it was shown that 7— law perfect fluid Bianchi type II models, 
are future asymptotic to the Collins-Stewart model [H] when | < 7 < y, consequently these 
models do not isotropise and the angle of tilt becomes negligible at late times. At the value 
7 = y the tilt destabilizes the Collins- Stewart model and there is an exchange of stability with 
the Hewitt model. In addition Hewitt solution has a four dimensional stable manifold and play 
the role of the future attractor for general Bianchi type II 7— law perfect fluid models. An 
interesting feature of the stability analysis given in (7j is that the equilibrium point 7 = y 

represents a second bifurcation between the equilibria (y, y) and (y,2) (the latter does not 
correspond to a self-similar model) and exhibits the property of the asymptotically extreme tilt 
for models where the state parameter 7 belongs to the interval fy, 2^. 

Concerning the remaining Bianchi models the only known solution is due to Rosquist and 
Jantzen (RJ) JOJ E] who have found a family of Bianchi VIo models in which the state 
parameter 7 lies in the interval (1.0411, 1.7169) (7 7^ 10/9). Although these models are known 
over a decade it seems that their importance were not fully understood. These models belong 
to the class n" = 0, they are rotating and consequently there are extra degrees of freedom 
adding to the difficulty of qualitatively analysing general Bianchi VIo models. Furthermore it 
was not known if this family represents the most general self-similar tilted perfect fluid type 
VIo solution. 

The overall goal of the present work is to fill this gap i.e. to show that the self-similar 
Bianchi VIq tilted perfect fluid solution given in ^T] is the most general but the state parameter 



7 G |, |J. In addition, we present in natural coordinates, the form of the self-similar metric, 
the Homothetic Vector Field (HVF) and the kinematical and dynamical quantities of the model 
and we discuss its physical implications in the asymptotic behavior of more general Bianchi VIo 
models. In the rest of the paper we follow the notations used in [6 . However for convenience 
we present briefly some of the basic results. 

In jH] the form of the transitively self-similar Bianchi VIq metric, the normalized tilted 
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fluid velocity and the HVF are given in terms of a set of arbitrary integration constants 
c a /3,v a ,a,b, D,ip (equations (3.29)-(3.31) of where ip is the homothetic factor. Setting 
a = — (pi + Vi)^ + 4^ — 6 and D = (2 — p-ijij) — b the frame components of the self-similar 
metric and the fluid velocity become: 



9a/3 



I Cll t 2 ( pi -V Cl2 t Pl+P2 - 2 C 13 t Pl 

Cl2t Pi+P2-2 C 2 2t 2{P2 - 1) C 23 t P2 
V c 13 t pi C 23 t P2 C 33 t 2 



A 1 = v 1 1? 1 - 1 ,A 2 = v 2 t p *- 1 ,A 3 = v 3 t (2) 
whereas the HVF assumes the form: 

H = iptd t + [(2 - p 2 )i; - b] d x + [-(pi + p 2 )V + 4^-6] yd y + fe<9 2 . (3) 
The FE and the Bianchi identities can be written: 

Rab - li^u a u b - 2 7 /lff a6 = (4) 

/i ;aM a + 7/20 = (5) 

( 7 -lfe fc + 7Ma = 0. (6) 

In contrast with the existence of two hypersurface orthogonal Killing Vectors (KVs) in type II 
models, self-similar Bianchi VIo models are necessarily rotational and admit only one hypersur- 
face orthogonal KV (X x or X 3 ) 8J. Therefore it is convenient to divide our analysis according 
to whether X^R^X^ = or X^R k ^ a X 3b \ = (since the KVs X 1; X 2 form an Abelian subgroup 
of the G 3 group of isometries the X 2 -case is similar). 
Case XfR k[a X lb] = 

We employ a new constant s which is defined by the relation: 

Pl +p 2 = 2(s + l). (7) 
Equation (j^J) implies that the state parameter 7 is related with s via: 

7 = . (8) 

[ 2s + 1 v ' 

From equation (£Q) (or equivalently the existence of the hypersurface orthogonal KV X x ) it 
follows that vi = and, in order to avoid the orthogonal case, we assume v 2 ,v 3 7^ 0. The 
resulting system consists of a set of 10 highly non-linear algebraic equations (@J) in 13 unknowns 
augmented by the system of equations (JHJ). 
Imposing the physical restrictions: 

fl > 0, T 2 > 0, a 2 > 0, uo 2 > 0, ifii a > (9) 
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we determine analytically the exact form of the self-similar metric and the fluid velocity. 1 

The rather lengthy and typical computations present no particular interest and we summa- 
rize the results below (without loss of generality we set C23 = — C\\ = 1): 

Cl2 = C13 = (10) 

C33 = {(p2-2)[p 2 + 2(s-l)]}- 1 (p 2 -s-l)[p 2 2 + 2p 2 (5s-2)- 
-4 (2s 2 + 3s - l)] x [3p 2 -2(2s + 1)] /[3pl - p 3 2 (s + 14) + 
+pj (23 - 10s 2 ) + 2p 2 (4s 3 + 18s 2 + 3s - 8) - 

-4(2s + 1) (2s 2 + 3s- l)] (11) 

c 22 = [3^2 - vl (s + 14) + p\ (23 - 10s 2 ) + 2p 2 (4s 3 + 18s 2 + 3s - 8) - 
-4(2s+l) (2s 2 + 3s-l)]/ 

{[32s (p 2 -s- l) 3 ] (p 2 - 2)- 1 \p 2 + 2 (s - I)]" 1 } (12) 

(2-»)fr + 2(«-l)] 
4(p 3 -s-l) 

= (j>2-s- 1) [3pj + 2p 2 (s - 4) - 4 (2s 2 - s - 1)] 

^ 3 (2-p 2 )[3p|-p|(s + 8)-p2(10s 2 + 2s-7) + 2(4s 3 + 8s 2 + s-l)]' 1 J 
The constant p 2 is related with the "state parameter" s according to: 



V4s 2 - 36s + 17 |3s - II - 42s 2 - 17 (s - 1) . . 

P2 = tz — 57 • 15 

17 — 3os 

The family of self-similar solutions ()10 )) -([15 )) has been given previously in ^T] using Hamiltonian 
methods and the full group of automorphisms of the Lie algebra of isometries. The advantage of 
the approach presented here is that we have used a natural choice for the coordinates, adapted 
to the canonical 1-forms {metric approach) hence we can study the resulting models directly. 
The kinematical and dynamical quantities of this family of self-similar models are: 

(2s + l)(2s + l-p 2 ) , . 

[i = - 2 (16) 

■ a- (P2-S- 1) [ P2 + 2 (s - 1)] (2s - l) 2 [3p 2 - 2 (2s + 1)] 

t 2 (p2-2s-l)[3p 2 -2(3s + l)] [U) 

r 2 = 3p 3 - p\ (s + 8) - p 2 (10s 2 + 2s - 7) + 2 (4s 3 + 8s 2 + s - 1) 

(p 2 -2s-l)[3p 2 -2(3s + l)] 1 } 



lr T\ie majority of the computations have been made using the algebraic computer package DERIVE. 
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a 2 = {3p 3 2 (56s 2 - 31s + 8) - p\ (56s 3 + 570s 2 - 249s + 64) - 

-Ap 2 (l40s 4 - 153s 3 - 161s 2 + 62s - 14) + (19) 
+4 (ll2s 5 + 90s 4 - 179s 3 - 40s 2 + 21s - 4)} x 

x {6t 2 [3pl - p 2 (12s + 5) + 12s 2 + 10s + 2] } _1 

= .fr + 2(.-l)]' 

2t 2 [3p 2 - 2 (3s + 1)] V ; 

The positivity of the above quantities is ensured provided that: 

It is easy to show that this solution belongs to the subclass of spatially homogeneous models 
satisfying the constraint = hence it represents a family of models which has a four 
dimensional stable manifold [12J. 

We also note that for the case where ^ there is a solution of the FE (@J) in which the 
constants c a p are all non-vanishing. For these models the state parameter 7 takes the values 
(l, |) U (|, 2). However none of these solutions is physically acceptable i.e. they do not satisfy 
the inequalities 0. 

L 3-Kfc[a^36] 



Case XkRuuXui = 



In this case v% = and the solution of the FE implies that P\= Pi = \ which by means of 
(J7J) and (JHJ) the state parameter 7 = |. The components of the self-similar metric and the fluid 
velocity are: 

C13 — C23 = 1, c 22 = Cn (22) 
9c 3 , + c 2 u (324c 2 2 - 387c 12 + 40) + c n c 12 (387c 12 - 80) - c\ 2 (9c 12 - 40) 
° 33 12 (c? 2 - c 2 n ) [ Cll (9c 12 - 1) + c 12 ] 1 dj 

.^-. 2 ^- 2[ % ( f 12 ; 1) + Cl2] (24) 

3 (en + C12) 

9c 3 x (54c 12 - 5) + c 2 u (810c 2 2 - 387c 12 + 40) + c n c 12 (441c 12 - 80) - c\ 2 (9c 12 - 40) = 0. (25) 

We note that this model also belongs to the subclass = and admits the hypersurface 
orthogonal KV X3. Furthermore it was recognized as an equilibrium point of the dynamical 
state space of the 7— law perfect fluid Bianchi VIo models [T5] . 

In conclusion we have shown that the general tilted Bianchi VI self-similar solution is given 
by (Q7J)-(Q3J) an d |Hf-(ll2P an d the state parameter satisfies | < 7 < |. It should be noted that 
in JT] the determination of the range of the state parameter 7 is based only on the positivity of 
the energy density /L Indeed it is straightforward to show that for 7 G (l.0411, |) U (|, 1.7169) 
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(7 7^ 10/9) the energy density is positive {fx > 0) but some (if not all) of the kinematical 
quantities (fTTjl - lj^Uj) are negative. 

The family of self-similar models found in the present paper could play a similar role (as 
in the case of Bianchi type II models) in the asymptotic behavior of general Bianchi type VIo 
models. In fact Barrow and Hervik ^3] have studied the (local) stability of tilted Bianchi 
models at late times (in the neighborhood of the non-tilted equilibrium point) and have shown 
that in type VIo, the non-tilted Collins solution [T3] (future attractor for nearby trajectories) 
is stable whenever 7 G f|, |). Therefore for 7 = | the tilt destabilises the Collins model i.e. 
there is a bifurcation which is associated with the stability change of the equilibrium points 
(Collins model and the present solution). In addition for 7 G (f> f ) ^ ne P resen t family of 
models is asymptotically tilted (t>3 7^ 0) and has a four dimensional stable manifold ^2] hence 
it is possible to play the role of the future attractor at least for the subclass of tilted Bianchi 
VIo models satisfying n° = (we recall that for general Bianchi VIo models i.e. when 7^ 
the dynamical state space is seven dimensional). 

In view of the above results and the similarities between types II and VIo one may conjecture 
that the subclass n™ = of Bianchi type VI models are possible to be asymptotically tilted 
at late times for 7 G ( |, 2^j. Furthermore a preliminary analysis on this class of tilted Bianchi 

VIo models indicate that whenever 7 G (§)2), the tilt is becoming extreme at late times [IS] . 

The author is grateful to Sigbj0rn Hervik for enlightening discussions and comments on an 
earlier draft of the manuscript and for bringing into the author's attention some of the results 
of his work. 
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